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POPULAR SUMMARY 

The Quasi Biennial Oscillation in the Earth’s upper atmosphere has an analog deep inside the Sun. As 
on Earth, the flow is east or west, it is at low latitude, and it reverses direction in a roughly periodic 
manner. The “period” in the solar case is 1.3 years. It was detected using solar oscillations similar to 
the way earthquakes are used to study the Earth’s interior. But its cause was not known. 

We showed that global oscillations (g-modes) can supply enough angular momentum to drive zonal 
flows with the observed reversal period. This required a calculation of wave dissipation rates inside 
each flow and in the turbulent layer that separates any two flows of opposite sign. Heat that this process 
leaves behind causes a thermal spike inside the Sun at the same depth. This may explain an anomaly in 
observed sound speed that has had no sure explanation. 
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ABSTRACT 

The reversing east-west flows centered near 0.675R (R=solar radius) cannot 
be explained by gravity waves generated by convection because' suitable wave- 
lengths do not penetrate deeply enough to drive the observed flow. An alternative 
is explored using g-modes, many of whose rotation rates have already been de- 
tected and are consistent with the known rotation of the nonconvecting interior. 
Angular momentum transfer to the fluid from a pair of asymptotic prograde and 
retrograde g- modes is derived in terms of their local dissipation rates and the 
flow velocity. This is used to show that a model flow with peak zonal velocities 
±200 m/s separated by a turbulent shear layer 3 Mm thick is consistent with 
the observations. The flow occupies the same thin layer where the sound speed 
is anomalous. Forcing this flow to replace itself every 1.3 years, as observed, 
requires power of about 0.3% of the solar luminosity. 

Subject headings: Fluid Dynamics — Sun: oscillations — Sun: interior 


1. INTRODUCTION 

Howe, et al. (2000) reported reversing zonal flows below the Sun’s convective zone (CZ) 
roughly periodic at 1.3 years. The ’’period” is not a constant (Antia & Basu 2000). More 
data, (Komm, et al. 2003) confirmed these motions and showed an anticorrelation between 
flows of comparable strength near 0.72R and 0.63R (R = solar radius). A correct model 
must have significant flow on either side of the midpoint, r m = 0.675R. 

Gravity waves generated in the CZ (by Reynolds stresses) and below it (by overshoot) 
were considered to drive the flow. Unfortunately, these waves fall far short of penetrating 
to r m and lack enough energy. Theoretical studies by Talon, Kumar, & Zahn (2002) and 
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Kim & MacGregor (2003) produce flows that are much too close to the base of the CZ 
compared to the observed flows. For example, when Kim & MacGregor (2003) computed 
yearly velocity reversals, they were < 0.005R, below the CZ base at 0.713R and rather slow 
<1 m/s, leaving unmodeled the great majority of observed flow volume and energy. Talon 
& Charbonnel (2003) also computed a very thin flow for a star similar to the Sun. When 
Talon, Kumar, & Zahn (2002) assumed a reasonable gravity wave power of 1.6 x 10~ 4 L 
(L = solar luminosity), they found a flow reversal period about 200 times too long. Since 
reversal time varies inversely as power input, their model would need almost 4% of the 
solar luminosity running downward as gravity wave flux in suitable wavelengths to give the 
observed period. This is very unlikely. 

A more penetrating mechanism with potentially more energy is explored here using 
mostly analytic methods. We propose that g-modes excited in the Sun’s core deposit 
angular momentum in the presence of these flows and drive them. Since the process 
deposits lots of heat, it can account for some or all of a conspicuous error in sound speed 
(§2), explaining why it is centered at the same place, r m , as the flow. In §3, momentum 
deposition from a pair of oppositely running g-modes is derived in terms of flow velocity 
and heat dissipation. This is applied in §4 to a model flow field whose diffusivities are 
calculated and cause the flow to drift downward by locally depleting the g-modes. A drift 
rate consistent with the observed 1.3 yr period can be achieved with reasonable numbers 
(<10 3 ) of small amplitude g-modes. 


2. SPIKES IN SOUND SPEED 


The fractional error Sc/c between predicted sound speed and helioseismic measurements 
exhibits two spikes (Basu et al. 1996; Kosovichev et al. 1997; Basu et al. 1997) approximated 
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on Figure la by halves of a cosine curve 

^ = Acos[k c (r - ri)] (1) 

centered at r* with amplitude A and wavenumber k c = 27r(3u>) -1 where w is full width at 
half height. The constants printed on the figure closely typify the last several figures of 
Basu et al. (2000). Both spikes ride on a broader discrepancy (wide arc) whose detailed 
shape is not needed here. 

Since c 2 is proportional to temperature T 0 divided by mean molecular mass /r, some 
have altered [i since the composition of solar models is fairly flexible at the percent level. 
Brun, Turck-Chieze, & Zahn (1999) and Elliott and Gough (1999) changed the helium 
mass fraction Y by about 0.01 and assumed mixing well below the CZ to make // and 5c/ c 
smaller. 


But T 0 must also be changed since reversing flows imply heat deposition. For an upper 
limit, let a spike be entirely due to a deviation T from model temperature To- Then T/T 0 
= 25c /c and the radiative diffusion equation 


dT cPT 
dt ~ T dr 2 


( 2 ) 


with equation (1) determine decay time T/T = (D r k 2 ) -1 . Figure lb gives D r for a standard 
solar model (Guenther, et al. 1992) and decay times are 7600 yr and 2.2 Myr respectively, 
for outer and inner spikes. The power to supply these radiative losses is roughly the excess 
thermal energy f dm(k b T)[('j - l)/r] _1 divided by the decay time, where k 6 is Boltzmann’s 
constant and dm is the mass element. Integration over the half-cosine on figure la and 
between latitudes ±20° occupied by the flow (presumed here to cause the spike) shows that 
about 0.05L must be deposited in the outer spike to maintain it in this example. 
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3. DRIVING A FLOW WITH g-MODES 

The g-modes (Cowling 1941; Christensen-Dalsgaard & Berthomieu 1991) are globally 
coherent oscillations trapped in the nonconvecting interior and evanescent above the base 
of the CZ. They drive zonal flows the same way as gravity waves but have a much greater 
source of energy. Wolff (2002, 1983) has already detected the unique sequence of solar 
g-mode rotation rates with high confidence (accidental probability, 5 x 10~ 4 ). Given 
this evidence that many angular harmonics are active, thousands of modes would not be 
surprising. The number of possible modes is i m ax(^max + 2) An for all angular states, 1 
and An radial harmonics for each angular state. One has 5,760 g-modes for 
£ max = 30 and An = 6. 

Modes concentrated toward the equator (\m\ > i/2) are of more interest since most 
flow energy is at low latitudes. Simple asymptotic properties will be utilized by assuming 
modes with n ^ l. The g-modes are excited in the nuclear burning core and transport 
net energy upwards to cover losses near and in the CZ. They deposit momentum in the 
same places, driving zonal flows and accentuating their shear (Ringot 1998). Flow will 
migrate downward toward the incident flux. Talon, Kumar, & Zahn (2002) and Talon 
& Charbonnel (2003) computed and plotted clear examples of this steepening shear and 
migration, which is upwards in their case because they use gravity waves carrying energy 
down. Of course, meteorologists have known all this for many years and call it the quasi 
biennial oscillation (Bretherton 1969), (+ several more REF’s). Mayr, et al. (2000) provide 
a recent model. 

In a uniformly rotating star and co-rotating coordinates, a linear g-mode with 
oscillation frequency a 0 has time-longitude dependence exp[i(m< 2 > — 0 - o t)] and zonal phase 
speed v P h = (< 7 o/m)r sin# with sin# set to 1 herein. A mode transports an upward flux of 
angular momentum Fj = (2m/ o^Fe where Fe is the upward energy flux (Zahn, Talon, & 



- 6 - 


Matias 1997). The doppler shifted oscillation frequency a = cr 0 (l-v/v p/l ) must be used 
instead of cro in any thin layer with zonal velocity v relative to the coordinates. Consider 
such a layer with effective horizontal area A in which g-mode energy is dissipated, causing 
a change A F E across the layer. Then the fluid gains energy at the rate, E — —AT# A and 
angular momentum at the rate, J = — AFjA. As a result 

dj _ 2m dE , . 

dt o dt' 

Now, as Kim & MacGregor (2003) and others have done, consider two modes that are 
identical except for the sign of m and dissipation rates that differ slightly due to doppler 
shifting and somewhat different amplitudes. Superscripts (+ and -) identify, respectively, 
the effect on the fluid of the prograde mode (m > 0) and retrograde mode (m < 0). The 
layer gains net momentum J 2 = J + + J _ , which by equation (3) and doppler shifting 
becomes 

= ^M[/3E 2 + E + - E-}( 1 - P 2 )- 1 (4) 

where ft = v/|n ph | and E 2 = E + + E~. For flows in this paper, /3 2 < 1. Since E 2 is 
always positive, its first term steepens shear by adding positive momentum to positive flows 
and negative momentum where v is negative. The remaining two terms cause downward 
migration due to ’’wave filtering” (§4). The single pair of modes discussed here stands for 
the effect of hundreds or thousands of such pairs. Modes with large |m| are most important 
for momentum transfer by eq. [4] and also because there can be many more of them. 


4. A MODEL FLOW 

The solid curve on Figure 2 models the flow as a radially damped sine wave at the 
instant its distance s from the Sun’s center is drifting lower past r = r TO . Since g-modes 
steepen shear as in papers cited above, we assume a fully turbulent layer develops with 
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thickness h = 3 Mm, bounded by velocities v s = ±120 m/s. Outside this ’’shear layer” peak 
velocities are \j = ±200 m/s achieved with the formula v = (355 m/s) sin(27rAx/0.04) 
exp(— |Ax|/0.015) with Ax = (r— ro)/R. Such a thin flow is hard to resolve. The gaussian 
G(r, ro) = exp(—|( Ax/0.032) 2 ) approximates current observational resolution, shown on 
Figure 5A of Howe, et al. (2000). Both have the same width, 0.075R, at half height. 
Observers can only see the convolution f drG(r,ro)v(r) divided by the area under G. This 
is plotted as a dashed curve, vertically magnified by 10 for clarity. It agrees reasonably with 
the Howe, et al. (2000) measurement, which has maximum power at 0.63R and 0.72R with 
its two peak velocities averaging 5 m/s. 

As the shear layer migrates down and weakens, an opposite shear will grow above it in 
an endless sequence of reversals. We presume from the observations that flows are strongest 
as a shear layer drifts past r m . Kinetic energy in the plotted flow is easily computed and 
has to be replaced about every 1.3 yr as turbulence in successive shear layers converts flows 
into heat. As a result the flow itself consumes a mean power < Ef > t — 0.0027L. 

Although tlae entire g-mode has virtually constant amplitude since the vast majority 
of its energy lies below the flow, its local amplitude is depleted (<10%) in traversing the 
shear layer or one of the unidirectional flow layers. We show this with mixing length 
approximations. For the shear layer, the overturning time = h/2v s = 3.5 hours and the 
diffusion constant D « 2v s h/3 = 2.4 x 10 8 m 2 /s. Then a g-mode with a 0 T ov < 1 would 
see its local temperature perturbation decay on a time scale = (Dk 2 ) -1 = 12 days if its 
<f> dependence dominates the effective wave number, k = m/r, and m = 30. For a vertical 
wavelength large compared to flow thickness, <8% of local oscillation energy is lost in an 
oscillation period < 1 day. Outside the shear layer, flow is vertically stable but it drives 
horizontal turbulence (Zahn, 1999) due to its limited latitudinal extent. One gets T m « 
(330 Mm)/v > 19 days, and D2 5 (330 Mm)v < 3.3 x 10 10 m 2 /s. In this case cjqt^ > 1 and 
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this factor lengthens the decay time to r d ~ a 0 T ov (Dk 2 ) 1 = (10 days)|v//v| 2 . 

Wave filtering causes downward migration of the flow. For brevity, let the two g-modes 
have the same amplitude and kinetic energy /mass, W / , where they are incident on the 
big negative flow (r = 0.655R). Due to their differing dissipation rates they are incident 
on the shear layer with energy/mass = (W f/Ha 0 ) f dr( 1 P)/Td. Integration is over 
the thickness H of the flow layer. Their difference drives the migration and is proportional 
to f drv 3 . Its value is W s + — W~ = 0.042 P where Pd is the oscillation period in days. 
When the shear layer on figure 2 moves a small distance to s— ds, a volume Ads with 
velocity — v s has been removed and a similar volume with +v s has been added. Fluid 
angular momentum is thus changing at the rate j = — 2v s rpAds/dt. With equation (4), 
the instantaneous migration rate of the shear layer is 

ds H (E+ -is-) 

dt a 0 pAv s r(l — f3 2 ) 

and (E + — E~)/pAh = (W+ — W~)/Td{ya. shear layer) to first order in f3. The numerical 
result is that ds/dt has the proper size to traverse the 0.06R thickness of the sound speed 
error in the required time ~1.3 years for an array of 1000 g-mode pairs with modest velocity 
amplitudes of 4 m/s below the CZ. Most of these modes would be undetectible at the solar 
surface. 


5. SUMMARY 

A plausible zonal flow field has been studied that agrees with flow observations by 
Howe, et al. (2000) and subsequent work. Only small amplitude g-modes are needed to 
drive the flow as long as at least ~ 10 2 mode pairs are active. This proposal conflicts with 
no other solar observations; indeed, it supplies a new explanation for the location of the 
spike in sound speed and lends support to earlier detections of g-modes. 
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The flow is not unique. A range of parameter space is available to fit flows of different 
strength and vertical wavelength. We are currently integrating more complete equations 
to constrict this flexibility and give better tests of such models. Although most people 
expected that, if the flow moved, it would move up towards an assumed gravity wave 
driver, Komm, et al. (2003) state that their work ” gives the visual impression that some 
fluctuations move ...[down]. ..into the lower convection zone and into the radiative interior”. 
Our g-mode driver would predict this. 
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Fig. 1. — a) Analytic fit to the fractional error 5c/ c, between sound speed from solar models 
from helioseismology . Narrow spikes whose amplitude, half width, and location are printed 
may be associated with deposition or extraction of heat by g-modes. The positive spike is 
centered at the same place as the east- west flows of Howe, et al. (2000) suggesting that flows 
and g-modes might be a partial cause. The negative spike at 0.18R is where Wolff (2002) 
inferred reflection of g-mode energy using their rotation rates, b) The radiative diffusion 
coefficient in a standard solar model also has a sharp peak just below the convection zone. 


Fig. 2. — The solid curve shows two oppositely directed zonal flows separated by a turbulent 
shear layer of thickness h = 3 Mm located a distance s from the Sun’s center. It is seen 
somewhere near peak strength as the system migrates deeper into the Sun and dies out. The 
flow pattern will be replaced from above by one of opposite sign about twice every 1.3 yr. 
When the model flow is viewed at today’s low resolution, the dashed curve results. It agrees 
with the observations in amplitude (5 m/s) and almost so in the location of its two extremes. 



0.004 



( s / 2 v w ) AjjAjsnjjia 



00 

o 


CD 

o 


" 3 - 

o 


CM 

o 


o 

o 


o/og 


Fractional Radius, r / R Fractional Radius, r / R 




Model Flow 



(s/lu) a ‘Ajpo|0A leuoz 



